Generating and Analyzing Constrained Dark Energy Equations of State and 

Systematics Functions 
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Some functions entering cosmological analysis, sucii as the dark energy equation of state or system- 
atic uncertainties, are unknown functions of redshift. To include them without assuming a particular 
form we derive an efficient method for generating realizations of all possible functions subject to 
certain bounds or physical conditions, e.g. w £ [—1, +1] as for quintessence. The method is optimal 
in the sense that it is both pure and complete in filling the allowed space of principal components. 
The technique is applied to propagation of systematic uncertainties in supernova population drift 
and dust corrections and calibration through to cosmology parameter estimation and bias in the 
magnitude-redshift Hubble diagram. We identify specific ranges of redshift and wavelength bands 
where the greatest improvements in supernova systematics due to population evolution and dust 
correction can be achieved. 



I. INTRODUCTION 

The nature of the dark energy accelerating the cosmic 
expansion is a major mystery of modern physics. The ef- 
fectively negative pressure giving rise to acceleration can 
be parametrized through the equation of state, or pres- 
sure to energy density, ratio of the dark energy. Obser- 
vational quantities such as the distance-redshift relation, 
Hubble expansion rate, or matter density perturbation 
growth (assuming general relativity) can then be derived 
in terms of the equation of state (EOS). However, little 
guidance exists from theory for the form of the EOS. 

One of the standard approaches is to adopt a well- 
tested, nearly unbiased functional form for the EOS. 
For example, the EOS as a function of scale factor, 
w{a) — wq + Wa{l — a), where wq, Wa are parameters 
to be fit, has been shown to be accurate at the 0.1% level 
in the observable distance for a wide array of dark energy 
models [l[. However, one may prefer to keep the EOS as 
free as possible. Values in bins of redshift, or some form 
of eigenmodes or principal components, do not impose 
assumptions on the form of w{a) (see, e.g., [1-0|)- 

Physics does bound the possible behaviors, though, not 
allowing full freedom in the bin values or principal com- 
ponent coefficients. One example involves a minimally 
coupled, canonical scalar field, where for all redshifts the 
condition must hold that w G [—1,4-1] for positive energy 
density. Principal components can be applied to many 
situations, such as the cosmic reionization fraction his- 
tory or the fraction of a source population in a particular 
subclass, where values can only lie in [0,1]. We gener- 
ically call such unknown, redshift dependent quantities 
"state functions" . Some approaches to such situations of 
"freedom under constraint" exist in the cosmology liter- 
ature, e.g. 0, i5i, iSi-ilOj. but here we concentrate on full 
and computationally efficient solutions. 

Going further into the motivation, we consider three 
reasons for imposing bounds: physicality, efficiency, and 
prior information. Some physical bounds are absolute. 



such as an ionization fraction ranging between and 1, 
while others are more relative, such as the dark energy 
equation of state ranging between —1 and -f 1 for a canon- 
ical, minimally coupled scalar field. In fact, there is a cer- 
tain amount of framework dependence in any analysis - 
the matter density cannot be less than zero, but the effec- 
tive matter density can appear less than zero when a uni- 
verse with a cosmological constant is interpreted in terms 
of a pure matter universe: this is precisely how the accel- 
eration of the universe was discovered. So while physical 
bounds are generally valid, results pushing up against the 
bounds should sound a note of caution; one might then 
loosen the bounds to check for consistent results. But 
starting with overly loose or unmotivated bounds has the 
price of computational inefficiency; in the vast majority 
of cases one would not scan over a space where the ion- 
ization fraction ranged from —5 to -1-5, say. Finally, the 
bounds may arise from prior information such as having 
measured a calibration offset to be less than some value 
(as we apply in Sec.|V|. There is little point in examining 
the effect of larger variations than allowed by this prior 
information. These rationales for bounds on the state 
function then translate directly into the principal com- 
ponent space. We emphasize that only the amplitude, 
not the freedom in the functional form, is being limited. 

When selecting physically valid principal component 
contributions the two main issues are those of purity 
- every set of values gives a valid state function - and 
completeness - every possible valid state function is rep- 
resented in the selection. In Sec. |ll] we discuss possible 
methods for generating principal component realizations 
of the EOS (or any other) function and assess their purity 
and completeness, especially when only a subset of modes 
is retained. We present a solution for the optimal - pure 
and complete - prescription in Sec. IIIIl along with an 
efficient mathematical shortcut and visualization for im- 
plementing it. We then turn to state functions represent- 
ing systematic uncertainties, whose evolution can cause 
incorrect cosmological conclusions. In Sec. IIVI we con- 
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sider supernova population fractions as the constrained 
state function and investigate the biases this can impose 
on cosmological parameters, and how to best constrain 
these with redshift specific observations. We discuss dust 
extinction corrections and their interaction with filter cal- 
ibration errors in Sec. [Vj and how to control these with 
wavelength specific measurements. The summary is pre- 
sented in Sec. |Vll 

II. REALIZATIONS OF THE EQUATION OF 
STATE 

We begin by phrasing the analysis in terms of the dark 
energy equation of state, although the results are gener- 
ally applicable to any state function. 

One possible goal for propagating an array of equa- 
tion of state functions to observational constraints is to 
place as little prior constraint on the functions as pos- 
sible, an admission of maximal ignorance in the hope 
that the observations impose form on chaos. A more re- 
strained approach is to treat the form of deviations from 
the basic function as free, perhaps representing unknown 
systematic uncertainties, though bounded in amplitude 
in some way. The most direct approach then is to de- 
scribe the deviations w{z) — Wb{z) by some value in each 
small redshift bin, equivalent to expanding in a top hat 
basis. 

This can be transformed into any other orthogonal ba- 
sis and we can hope that a principal component analysis 
lets us compress the information in some way, such that 
a small, tractable number of modes gives a simplified, 
though still somewhat diverse, functional form. We can 
write 

w{z) - Wb{z) ^^aiCiiz) , (1) 

i 

where we refer to w{z) — Wb{z) as the state function, 
ei{z) as the modes or principal components, and at as 
the mode coefficients. Note that the state function is re- 
ally the deviation from some baseline, and can represent 
the dark energy equation of state or the cosmic ionization 
fraction, supernova subclass population fraction, etc. In 
the top hat basis, ei{z) would simply be 1 within the ap- 
propriate redshift bin and outside, and a; would simply 
be Wi, the value of the state function within the bin. 

The state function may not be allowed to have ar- 
bitrary excursions, but can be constrained by physical 
or theoretical expectations to lie within some bounds. 
These bounds could be elementary, such as the ioniza- 
tion fraction must lie between and 1, or more physical, 
such as the equation of state for a minimally coupled, 
canonical scalar field must possess w{z) G [— We 
define the bounding function, or envelope, by 

W-{z) <w{z) ~Wb{z) <W+{z) . (2) 

Given real data, the results should localize within the 
bounds. One might be tempted to loosen the bounds and 



allow the data to lead to the proper area of parameter 
space. However, the data does not always have the re- 
quired leverage to make this a successful approach. For 
example, if the equation of state rapidly oscillated be- 
tween — 10 and -1-10, this could not be detected in the 
distance measurements (or such an oscillation in ioniza- 
tion fraction, even to unphysical negative values, might 
not be seen in cosmic microwave background polarization 
measurements) but one has spent a lot of effort calculat- 
ing over an enlarged range. Furthermore, when dealing 
with systematics, unknown by definition, or projected 
future measurements, if one does not bound the ampli- 
tudes then no real information can be obtained from the 
results. Thus, one has to balance reasonable, physical 
bounds and the computational efficiency with the desire 
not to restrict the input. We take W-, W+ to be defined 
with this in mind. The effects on the principal com- 
ponents of increasing the envelope are simply given by 
scaling W-, W+ in the formulas derived. 

The question then becomes how to best incorporate 
these bounds in "configuration" (e.g. redshift) space into 
the coefficients of the principal components (PC) in mode 
space. For example, to generate realizations of state func- 
tions that are viable according to the bounds imposed, we 
must know how to properly sample the PC coefficients. 

Several methods can be attempted, but must be as- 
sessed for their (computational) efficiency, purity, and 
completeness. The obvious, and least efficient method 
is simply to try values of the coefficients {ai, . . . , ajsi} 
and see if the functions wlz) obey the bounds at each 
redshift. In practice one must truncate the number of 
PCs at a finite number, choose a finite range for each a^, 
and with a number of grid points R sampling the coeffi- 
cient range evaluate i?^ functions to test whether they 
lie within the bounds. For a grid of 20 points and 10 
PCs, this requires 20^" « 10" evaluations. We call this 
the scanning strategy. It would be pure and complete, 
but is not efficient. 

A second approach is to ask that each PC contribution 
to the state function obey the bounds individually. This 
mode-by-mode strategy has been implemented in f^, Q 
for example. Projecting a given PC against the state 
function yields the coefficient: 

tti = TV J dz ei{z) [w{z) - Wb{z)] , (3) 

where J\f — 1/ J dzef{z) is the normalization factor. In- 
corporating the bounds on the state function and break- 
ing the integration region into those redshifts where (z) 
is positive and those where it is negative, one obtains the 
bounds 

a- <a,< a+ , (4) 

where 

a± = (A/-/2) J dz { [W+{z) + W^{z)] e,(z) 

± [W+{z)-W^{z)]\e,iz)\}.{5) 
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The main problem with this approach is that the 
modes are treated independently. So if one saturates the 
bounds on each coefHcient, say, then the generated state 
function may actually lie outside the envelope. Thus, this 
method is complete but not pure. 

One way to incorporate all the mode information is to 
consider the integral of the square of the state function 
0. Then 

J dz [wiz) ~ wb{z)]^ = ■ (6) 

i 

Imposing the bounds on the state function then delivers 
the constraint 

J2o^^ < (1/A^) / dzmax{Wl{z),Wl{z)}, (7) 



III. A PURE AND COMPLETE 
PRESCRIPTION 

The physical bounds on the state function are imposed 
in the redshift space but we need to translate these into 
PC coefficient space if we want to generate principal com- 
ponent analysis (PCA) realizations of the state function. 
The problem is that a principal component contributes 
to the state function over the whole redshift interval con- 
sidered. Effectively, PCA mixes the values Wi from all 
redshift bins in the bin basis. Therefore what the chosen 
bound corresponds to ~ a value for the state function at 
a particular redshift 2: - is not localized in coefficient pa- 
rameter space but is described by a linear combination 
of many modes weighted by their respective coefficients. 



A. Hypersurface Picture 



where the maximum is to be evaluated for each redshift. 
(This generalizes the expression in Q to when the enve- 
lope is not redshift independent.) We call this the inte- 
grated method and it defines a sphere in the mode coeffi- 
cient space. This approach guarantees completeness but 
not purity, i.e. every viable state function can be gener- 
ated with this set of coefficients, but nonviable ones can 
be as well. By itself it lacks a specific prescription for 
implementing the selection of a^'s. 

An alternate approach is the "global" method, where 
coefficients are chosen based on the previous coefficient 
values. For example, choose the coefficient ai based on 
the envelope constraint as if this were the only mode, 
giving 

ei[z) ei(z) 

where the top (bottom) sign holds for ei(z) > (< 0). 
This is applied for all redshifts under consideration and 
the tightest constraints obtained define the range of ai. 
Once an ai is sampled within the allowed range, one 
obtains similar bounds on a2 using, e.g., 0262(2;) < 
W+{z) — ai ei{z), and so on. This global method is sim- 
pler, not involving any integrals, although it still involves 
scanning over choices for the coefficients within their al- 
lowed range. However, a value for the coefficient at that 
has been rejected because it lies outside the bounds of 
Eq. ([5]) or similar may actually be valid because another 
PC counteracts its contribution and pulls the state func- 
tion back within the envelope. Thus the method is pure, 
i.e. all generated state functions will be viable, but not 
complete. 

Thus we have generating methods that are pure and 
complete but inefficient (scanning method), complete 
(mode- by- mode method) , and pure (global method) , but 
no obviously optimal method. We address this lack in the 
next section, and show how all the methods are related. 



However, by making use of the properties of linear 
transformations, and a particularly clear geometric pic- 
ture, we can implement an exact, fast method for the 
translation. Consider the envelope on a single Wi. This 
gives a range, or line segment, along the Wi axis. Combin- 
ing the envelopes for all redshift bins, i.e. Wi parameters, 
defines a hypersurface in an N dimensional space, where 
N is the number of redshift bins. If the bin bounds do 
not depend on values Wi in other bins, i.e. each bin is in- 
dependent (recall the original motivation was to consider 
state function behaviors without assuming a functional 
form), then the surface is a hyperrectangle. 

The corners of the hyperrectangle are defined by the 
values W±{zi) of the envelope. These 2^ vertices contain 
all information on the boundary between the permitted, 
i.e. viable, instances of state functions w{z) — Wb{z) and 
the disallowed or unviable ones. That is, the boundary 
defines the pure and complete set. 

We defined the vertices as sets of Wi coordinates but 
now let us consider the hypersurface in the PC coefficient 
space of ai coordinates. Because the PCA is a (nor- 
malized) linear transformation of the redshift bin values, 
the hypersurface is merely rotated, not distorted or ex- 
panded. If we are interested in a subset of M modes, 
smaller than the maximum number N (there cannot be 
more modes than the original bins used to define the 
PCs), then this corresponds simply to a projection of the 
hypersurface onto the subspace of the M PC coefficients. 
We illustrate the case of a 3 dimensional hyperrectangle 
projected onto 2 PC coefficients in Figured] 

The boundary of the pure and complete set of PC coef- 
ficients is defined by connecting the outermost projected 
vertices ensuring the boundary remains convex. This fol- 
lows from the linearity of the transform: the hypersurface 
is convex and so the projection must then itself be convex. 
(Note the projected figure is not in general rectangular.) 
The projection can be computed quite quickly through 
the use of matrix algebra (see the Appendix). Thus a 
pure and complete set of PC coefficients for viable, and 
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FIG. 1: A cube in bin space, corresponding to redshift in- 
dependent limits on w{z) — Wb{z) in 3 redshift bins (i.e. the 
simplest nontrivial case of an TV-dimensional hyperrectangle) , 
is projected onto the parameter space of the coefficients of the 
first two principal components. The small blue dots indicate 
the projection of the cube vertices. The permitted area in the 
ai-a2 plane is indicated by the light grey shading. Note that 
it does not form a rectangle. 



only viable, state functions can be generated efficiently. 

If we do a complete projection over all N dimensions 
except one, say aj, we obtain the absolute minimum 
and maximum bounds on aj. These bounds are equiv- 
alent to those found with the mode-by-mode method, 
Eq. Q. Furthermore, since distances are conserved un- 
der the linear transformation, the maximum distance in 
bin space, i.e. the longest diagonal of the hyperrectangle, 
must also be the maximum distance in PC coefficient 
space. A (hyper)sphere with this diameter circumscribes 
the allowed set and corresponds to the sphere of the in- 
tegrated method, Eq. ([7|). From the fact that the hyper- 
sphere circumscribes the hyperrectangle, it is clear that 
this method generates a complete, but not pure, set. 

The degree of impurity or incompleteness for the meth- 
ods can be tied to the ratios of areas (or hypervolumes) 
between the geometric figure defined by the methods and 
the true hyperrectangle. If the PC coefficients are highly 
independent of each other (of course the mode vectors 
themselves are orthogonal) , then we expect the mode-by- 
mode method, where we ignored the effect of aj on ai, 
to be a good approximation, i.e. nearly pure and com- 
plete. Taking zero correlation between coefficients de- 
fines a rectangle in the ai-aj plane, for any a^, aj. As 
the coefficient parameters become more correlated, the 



mode-by-mode method should become less efficient at 
finding only the viable state functions, i.e. less pure. Ge- 
ometrically, the filling factor of the true hyperrectangle 
projection will decrease. 

Figure [5] illustrates this relation between correlation 
and filling factor. The top panel shows the projection 
onto the space spanned by the coefficients of PC modes 
1 and 17. Since ei and en have their main weights at 
very different redshifts, the coefficients ai and an are 
substantially uncorrelated, and indeed the filling factor 
is high (but not perfect). The bottom panel displays the 
equivalent projection for modes 1 and 2. Here the over- 
lap of the modes in redshift is greater and so the coeffi- 
cients are more correlated; the filling factor is noticeably 
decreased. Therefore the mode-by-mode method is not 
efficient when considering the dominant modes. 

By contrast, the exact hyperrectangle projection 
method is highly efficient. For 10 modes, say, there are 
2^" w 10'^ vertices to evaluate (the projection takes neg- 
ligible computational time using the method in the Ap- 
pendix). Contrast this with the previous 10^'^ evaluations 
needed for direct scanning. 

If we were to increase the number of redshift bins (i.e. 
bin modes, holding the redshift range constant), this al- 
lows for more and more PC modes. However, since most 
of these additional modes would be less and less corre- 
lated with a given mode, we effectively have a conver- 
gence in the behavior of the parameters, i.e. the projected 
boundary in a given ai-aj plane. 

In summary, we have presented an efficient, pure, and 
complete method of obtaining the boundary defining the 
set of viable state functions. The relation to previous 
(not simultaneously pure and complete) methods is illus- 
trated in Figs. The outer rectangle gives the pre- 
scription of the mode-by-mode approach; the thick inte- 
rior polygon shows the exact solution using the projec- 
tion of the hyperrectangle; and the light shaded interior 
non-rectangle illustrates the global method, representing 
a cut through the hyperrectangle at ai^2 — 0. Since 
the exact solution can be generated efficiently there is no 
need to use the over- (global) or under- (mode-by-mode) 
approximation. 

While we have solved the problem of obtaining effi- 
ciently the constraint on the region of principal compo- 
nent space that is viable given some bounds on the state 
function, we have to ask whether this is really the best 
path for analyzing the effect of various state functions. 
To scan over all viable PCs we would select from the PC 
coefficients within the allowed region. If the probability 
of the state function in the bin basis was uniform within 
the bounds, then because of the linearity of the principal 
component transformation the interior volume in PC co- 
efficient space can also be uniformly sampled. However, 
in general we would have some correlation 

(aiaj) = eip Cjq {wpWq) , (9) 

where angle brackets denote the ensemble average and p, 
q are redshift bin indices (implicitly summed) while i, j 




FIG. 2: Projections onto PGA coefficient planes are sliown for 
tlie 17 mode case. The liglit, red vertical and horizontal lines 
indicate the mode-by-mode bounds of Eq. Q. The dots give 
the projection of the 2^^^ vertices of the hyperrectangle, and 
the thick, green polygon gives the outer boundary, defining 
the interior region of pure and complete viable state functions. 
[Top panel] The ai-an plane in parameter space has low 
correlation between these modes, so the filling factor of the 
approximate bounds is high. [Bottom panel] The ai-a2 plane 
has strong correlations because the PGs overlap substantially 
in redshift, so the filling factor is low. 



are component mode indices. 
Writing this in matrix notation, 

A = EWE'^ . 



(10) 



where W is the correlation of the state function (e.g. 
equation of state values in redshift bins) and A is the 
correlation that then must be imposed on the selection 
of PC coefhcients. Note that when the bounds on the 
state function are redshift dependent - as when some 
data constraint knowledge is incorporated - then even a 
diagonal W does not lead to a diagonal A. 




Mode-by-mode 
True permitted 
Global method 



FIG. 3: Three different approaches to obtaining the principal 
component coefficients for a constrained state function are 
illustrated, here for the first and second coefficients and the 
simplest case with only 2 modes. The true viable region is 
the thick black rectangle in wi-W2. The shaded, green interior 
polygon shows the pure but incomplete global approximation, 
while the exterior thin, red rectangle shows the complete but 
impure mode-by-mode approximation. We do not show the 
circle circumscribing the outer rectangle that corresponds to 
the complete but impure integrated approximation. 




FIG. 4: As Fig. [3] but for 3 binned w modes projected onto 
the plane of the first two PG coefficients. Here we suppress 
the axes. The light, red vertical and horizontal lines show 
the mode-by-mode limits; the thick, black polygon defines 
the exact, pure and complete function space, and the green, 
shaded region is the global approximation. 
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B. Restricting Modes 

If we keep only M modes in PC space then, because the 
ai axes are not in general aligned with the bin basis Wi 
axes, the PCs will still span the redshift range but will not 
be able to describe the full range of w{z) behaviors within 
the true bounds. That is, we diminish the completeness 
if we restrict the number of PC modes. Note this can be 
treated in the hyperrectangle picture as slices through 
the TV-rectangle at fixed values of the neglected N — M 
parameters (see the Appendix for more details). Figure[5] 
shows an example of the diminished state function space 
accessed when limiting to 4 modes (out of 17). Moreover, 
the impurity of the mode-by-mode method becomes more 
severe, with Fig. [5] showing that using less than the full 
number of modes can yield up to 70% of the generated 
forms of the state function being spurious, i.e. ones that 
invalidly exceed the bounds, under the restriction to the 
first M modes. 



S 0.0 




FIG. 5: Restricting the number of PC modes kept in the 
analysis will lose completeness in the state function forms 
allowed. The light, green dots show the projection of the 
full set of hyperrectangle corners, with the light, red outer 
polygon giving the pure and complete bounds in ai-ai space. 
The black dots and black polygon show the case when only the 
first 4 modes are kept. The region in between the polygons 
represents viable, but lost state functions. 



In the end, then, because of the coefficient correlations 
and the completeness issues, little advantage accrues in 
fact to the use of PCA for the scanning over functional 
forms. It is more efficient (and innately pure and com- 
plete) simply to carry out the analysis in the original 
state function space where the constraints originated. A 
standard redshift bin basis allows the freedom needed to 
model the form of the state function, and the constraints 
can be imposed naturally without complicating the gen- 
eration of realizations. In the next sections we demon- 
strate the real world application of the bin basis state 
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FIG. 6: When not all PC modes are kept then the purity of the 
mode-by-mode method decreases. The curve shows the range 
purity - the product of the bounds on each PC coefficient for 
the exact method when only M modes are retained vs. when 
all modes are. This is also equivalent to the ratio of the M 
mode range area to the mode-by-mode method's area. 



functions to problems involving calculating the effects on 
cosmology results when confronted with unknown sys- 
tematics functions. 



IV. SYSTEMATICS: SUPERNOVA 
POPULATION DRIFT 

The use of constrained functions, and their impact on 
parameter estimation or the science results, enters into 
myriad areas of cosmology. This is a particularly im- 
portant issue for systematic uncertainties, where we do 
not know the form of the residual error function. We 
therefore consider the example of the population frac- 
tion of a certain type of source as a key element of the 
cosmology calculation and take as the state function the 
uncertainty in our knowledge of it. By definition the 
function is constrained to take values in the range [0, 1]. 
If the source is a standardized distance indicator such as 
Type la supernovae (SN) and we posit that the popula- 
tions represent subclasses with slightly different intrinsic 
magnitudes, then any variation with redshift in the pop- 
ulation fractions will appear as magnitude evolution and, 
if unrecognized, bias the cosmological parameter estima- 
tion. This is known as population drift (for theoretical 
discussion and observational limits see |lll4la |). 

In the effects of population drift as a bias or in- 
creased dispersion (if adding fit parameters) on cosmol- 
ogy were investigated for a class of state functions de- 
pending as a power law in redshift (also see [13 )■ Here 
we can analyze every form of population drift and investi- 
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gate which are the most dangerous. In addition we refine 
our quantification of the cosmology bias and explore in 
what redshift ranges the population drift systematic is 
most biasing. 

Population drift as a systematic relies on two elements: 
an actual difference in intrinsic magnitudes between the 
subclasses and a redshift dependence in the difference. 
A mere constant difference is absorbed into the absolute 
magnitude nuisance parameter A4. For simplicity, we 
illustrate the basic results for a two population model, 
where the SN have a fraction fi{z) with intrinsic mag- 
nitude Ml and a fraction /2(z) with intrinsic magnitude 
M2. We can consider /i as representing all the popula- 
tions we recognize and /2 = 1 — /i as an aggregate of 
those unrecognized. Then the unrecognized systematics 
appears as a magnitude evolution 

Am(z) = AAf [/(z) - /(O)] , (11) 

where AM — M2 — Mi and / = /2 is now our state func- 
tion. The constraint on the state function, by definition 
of the population fraction, is / G [0, 1]. 

Propagating this systematic through to the cosmology 
parameters is straightforward. For a parameter set p — 
{pi}, the bias is (see, e.g., fl8| ) 

Sp^{U^C-^U)-^U^C-^ AO, (12) 

where O is the observable, U — dO/dp, C the error co- 
variance matrix for the observables and AO is a system- 
atic offset in observable O. The term in parentheses is 
simply the Fisher matrix, and so its inverse is the pa- 
rameter covariance matrix. In the case we are currently 
considering the observables O are SN magnitudes at var- 
ious redshifts and AO is the magnitude offset of Eq. pT|) . 
For a diagonal error covariance matrix the equation takes 
a simpler form 

^ 1 flm, 

Sp. = AMiF-%Y.n^,)-^^^. (13) 

Here F is the Fisher matrix, say 4x4 with respect to A4, 
f^m, Wo, and Wa, where flm is the present matter density 
in units of the critical density. There are N data points, 
each with an associated redshift Zk- 

We can now explore the effect of any form for f{z), 
subject only to the constraint < f{z) < 1. We do not 
need to assume a functional form for /, rather we want 
to allow it complete freedom under constraint. As we 
saw in the previous section, no real advantage accrues 
to a principal component analysis - all the information 
exists and is more accessible using a redshift bin basis. 
In fact, PC A when keeping only a more limited number 
of modes loses information and the bin basis allows for 
greater efficiency in scanning the allowed state function 
parameter space. 

Before we calculate the bias we examine in more de- 
tail how to assess it quantitatively. The bias Spi itself is 



only informative together with the cosmological param- 
eter uncertainties. If the estimated uncertainty on the 
parameters is large, then the relative effect of a particu- 
lar bias is lower, meaning that the (mis)estimated model 
is still within some acceptable confidence level contour. 
For each parameter pi , Q employed the risk statistic \^ 

Risk(p,) = + Sp^ (14) 

as a measure of the influence of the bias. However, the 
overall cosmology is biased by the vector dp. One could 
imagine that each parameter bias relative to the disper- 
sion is small, but in a direction such that dp is oriented 
along the thin part (minor axis) of the confidence level 
contour; then a small shift could actually be a large bias 
relative to the contour, i.e. in terms of the Ax^. Follow- 
ing [20!] (cf. I^H) therefore, we use as our bias statistic 

Ax' =5pFM<5p^, (15) 

where dp is the vector of parameter biases we consider 
and F^*") is the reduced Fisher matrix, marginalized over 
all parameters except those in whose biases we are inter- 
ested. For example, if we consider biases in the WQ-Wa 
contour, then F^'") is the inverse of the 2x2 submatrix 
of the covariance matrix containing wq and Wa- In the 
case of a single parameter. Ax' = {5p/apY and Risk 

= apVl + Ax'. 

We can now scan over all possible population drifts 
f{z) — /(O) and evaluate the bias effects on the cosmo- 
logical parameters. We write / in the redshift bin ba- 
sis, initially with 17 bins uniform between z = Q — 1.7. 
For the Fisher matrix we take simulated data based on 
the SNAP SN redshift and error distribution plus 
a Planck-inspired constraint on the reduced distance to 
CMB last scattering of 0.2%. The fiducial cosmology is 
ACDM with matter density 17„ = 0.28. 

Table U describes the population evolution functions 
computed to deliver the maximum bias in Ax'. The re- 
sults have a very simple form: a single or double sharp 
transition in redshift. This can be understood through 
analyzing Eq. (jl3p . The bias is a linear transforma- 
tion of /, hence for a maximum bias / is driven to the 
extreme value that complements the sign of the term 
Pi = {F~^)ij dnik/dpj, for each Zk- To maximize Spi, 
when f3i > then / should be 1, while when /3i < then 
/ should be 0. This will give coherent addition of the 
terms in the sum and so deliver the largest Spi. Thus / 
should simply be a series of tophats over those redshifts 
where I3i(zk) is positive. For the parameter flm, /3r2„, 
crosses once through so SQm is maximized by a pop- 
ulation function that has a single transition, at z « 0.7. 
Thus the most potent evolution function / - the one 
having the strongest consequence for cosmology estima- 
tion - has a step appearing at z = 0.7 and extending to 
the maximum redshift. For wq or Wa, the respective /3's 
cross twice through so / forms a tophat extending from 
z « 0.2 (respectively 0.1) to z « 1.0. 
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Parameters ztrana max Axam=o.oi AM (la) 





0.7 


0.92 


0.010 


■Wo 


0.2, 1.0 


1.35 


0.0086 


Wa 


0.1, 1.0 


1.35 


0.0086 


Wo, Wa 


0.1, 1.0 


1.36 


0.013 


Qm,W0, Wa 


0.1, 0.9 


1.39 


0.016 



TABLE I: For each set of parameters we consider the form 
of population evolution that maximizes the cosmology bias 
in terms of Ax^. The ■St rans column eives the redshift of the 
maximizing step function in /, dehvering a maximum bias 
Ax^ scaled to the case where AM = 0.01, shown in the next 
column. Note Ax^ will scale as AM^. The last column shows 
the value of AM that will shift the derived cosmology by la 
from the true cosmology. 



This has a number of crucial implications. First, since 
the state function with the maximal effect arises from 
a sharp transition, we see that it was prescient to use 
the bin basis for the population function. Had we trans- 
formed to principal component space, or some smooth 
orthogonal basis such as the Chebyshev polynomials con- 
sidered by [i,[l^, then we would have had difficulty ap- 
proximating the true solution with a finite number of 
modes. Second, the sensitivity to population evolution 
at specific redshifts guides the survey design to obtain 
especially detailed measurements at these redshifts. The 
results indicate that as observations make the transition 
from local (2 < 0.1) SN to low redshift {z > 0.2) SN, they 
must comprehensively collect and study the SN proper- 
ties so as to ensure a firm like-to-like comparison and not 
allow for unrecognized populations. Similarly, the tran- 
sition from z < 1 to z > 1 SN is key, so a transition from 
ground-based observing to the space-based observing ne- 
cessitated at z > 1 could be problematic. A homoge- 
neous survey extending across this transition would have 
far better control over the systematic uncertainty. 

These results hold as well when considering simultane- 
ously bias in multiple parameters in the Ax^ formalism. 
The population / enters Eq. quadratically and so 
the bias is still maximized by the extreme values of /, 
i.e. top hats in redshift. The transition locations do not 
shift appreciably when considering bias in the two pa- 
rameter space of WQ-Wa nor the three parameters fim. 
Wo 7 Wa simultaneously. Furthermore, the transition loca- 
tions are robust to changing the step functions to more 
gradual slopes. 

The last column of TableUshows the magnitude of AAI 
that in the worst case of population evolution causes a Icr 
misestimation of the cosmology. This is where the scal- 
ing of the state function bounds enters: the shape, i.e. 
redshift dependence of the population function, is unaf- 
fected by amplitude of the bounds, but the absolute level 
is determined by the bounds. If we consider twice as large 
values for AM (or if we were to unphysically allow / to 
range from to 2), then Ax^ just scales with AM^. If 
we want to be sure that population drift cannot cause a 




0.000 0.005 0.010 0.015 0.020 0.025 
AM 



FIG. 7: For the maximally biasing population drift the cosmo- 
logical parameter set {flm,Wo,Wa} is biased by la for a pop- 
ulation magnitude difference AM = 0.016 mag. The curve 
shows by how many a the best fit cosmology is biased as a 
function of AM. 



> Icr shift in the equation of state parameters, we need 
to be able to recognize SN subclasses differing by 0.0086 
mag or more. 

Figure [7] shows the relation between the maximum 
number of standard deviations a by which the cosmology 
is distorted, as a function of difference in absolute mag- 
nitudes AM between the populations, for the full set of 
cosmological parameters. While Ax^ scales as AAf^, the 
number of a this bias corresponds to only scales as AM 
in the one parameter case. We see that for three param- 
eters the a remains nearly linear for large AM but does 
not improve as rapidly for AM < 0.015. 

Beyond the maximum Ax^, we can investigate other 
properties of the biasing. For example, we can explore 
further the direction of the systematic shift caused in the 
cosmology parameters, the relation of the forms of the 
population drift, i.e. the number of steps or oscillations, 
to the bias, and the overall statistics of the biasing. 

We begin with the effect on the equation of state esti- 
mation caused by the systematic error. Figure [8] shows 
the specific form of bias induced in the equation of state 
by the 10 worst case population drifts. The worst biases 
all distort the cosmology in the same way: making a cos- 
mological constant look like a rapidly varying equation of 
state. Indeed, this is characteristic not just of population 
drift but of any sharp transition in the SN magnitudes, 
such as from patching together two redshift samples with 
an unrecognized offset (local to low redshift samples, or 
ground-based to space-based). This points up the need 
for tight crosscalibration, and ideally a continuous, ho- 
mogeneous data set, as well as the need for caution in in- 
terpreting an apparent behavior of the equation of state 
crossing w = —1: exactly what is expected from such a 
systematic. 
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FIG. 8: The biases induced for AM = 0.02 by the worst pop- 
ulation evolution functions - those that induce the 10 largest 
Ax^ - are plotted in the wo-Wa plane as x's. Note how they 
congregate at the extreme end of the major axis of the confi- 
dence contour. The inner ellipse indicates the 39% confidence 
level contour (so a{wo), a'{'Wa) are given by direct projection 
to the respective axis), while the outer ellipse shows the 95% 
cl joint likelihood contour. 



The influence of the forms of the population evolution 
function on the bias generated in the cosmology parame- 
ters can be investigated through looking at the statistics 
of the Ax^ distribution. For example, while the maxi- 
mum bias is generated from a population function with 
one or two steps at sensitive redshifts, we expect a large 
number of transitions to have relatively little effect since 
such an oscillatory behavior does not resemble the effect 
of a cosmological parameter. The sum of the terms in 
Eq. (IT3)) effectively cancels out. Thus certain types of 
systematics are fairly benign, such as quasi-periodic k- 
correction errors O [23| • Figure IHl shows the range of 
Ax^ generated as a function of the number of transitions 
in / between redshift bins. As expected, as the number 
of transitions gets large, the bias decreases. Similarly, 
when the step amplitude is small then the bias is negli- 
gible so the Ax^ distribution ranges between and the 
maximum for each number of transitions. 

The location in redshift of the features in the popula- 
tion function also are important. As we saw in Table HI 
z « 0.1 and z « 1 were key regions for sensitivity to bias. 
In Fig. [TU] we plot the redshift locations giving not just 
maximum bias, but greater than a certain percentage of 
maximum bias (still keeping full steps, i.e. / = or 1). 
We see that down to 50% of the maximum possible bias 
the culprits are still population evolution around these 
sensitive redshifts. This suggests that surveys designed to 
recognize population subclasses through especially com- 
prehensive measurements around these redshifts can re- 
move the top half of possible cosmology bias, improving 
the systematics by a factor two. 



FIG. 9: The range of possible biases Ax^ in the cosmology 
likelihood contour for flrn-wo-Wa is plotted vs. the number of 
transitions (changes in value between the 17 redshift bins) in 
the population evolution. 

If we consider every possible form of the population 
function, randomly scanning over the number of steps 
and locations of transitions, then most of these will have 
little effect on the cosmology. The mean bias, or (Ax^) 
will therefore be small. For example, the mean is only 
0.08 for the wo-Wa contour. This is simply due to com- 
binatorics: there are many more ways of having, say, 8 
steps over 17 bins than 1 step - some 24000 times more 
possibilities - and multistep functions will have little im- 
pact on the cosmology. Furthermore the random loca- 
tion of the steps will also dilute the mean bias. However, 
random population evolution is not the issue; for survey 
design we have to consider the worst case scenario, i.e. 
which systematics can give the most egregious misesti- 
mation of the cosmology results, and how to control this. 
The results indicate that experiments should be guided 
by requirements to recognize subtypes with magnitude 
differences down to ^ 0.01 mag (Fig. [7]), and with par- 
ticularly comprehensive measurements around z w 0.1 
and z « 1 (Fig.Iini). 

V. SYSTEMATICS: DUST CORRECTION AND 
CALIBRATION 

The analysis of the constrained state function in terms 
of population drift was particularly straightforward be- 
cause of the linear relation between the function /(z) 
and the observable m{z). To illustrate a more compli- 
cated application we consider the systematic uncertainty 
due to dust extinction correction in supernova distances. 
This is currently one of the dominant systematics [25l - [29| 
and uses measurements in multiple wavelength bands, or 
filters, to correct for the dust effects. However, if the dif- 
ferent filters have some uncertainty in their calibrations 
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FIG. 10: The histograms show the redshift locations of the 
most sensitive steps in the population function / for various 
cuts in Ax^. In order from highest to lowest peak the cuts 
are 90%, 75%, 50%, 25% of the maximum possible Ax^. The 
top panel considers the Ax^ bias for the cosmology likelihood 
contour in the Qrn-wo-Wa space while the bottom panel is for 
the WQ-Wa space. Note that the most sensitive redshifts are 
robust to lower levels of bias, down to 50% of maximum, and 
to the parameter space considered. 



then this propagates through to the relative fluxes or col- 
ors and then to the dust correction [s^l- We use a simple, 
two band version of this as an illustration of a nonlinear, 
constrained systematic. 

We take the systematic to arise from zeropoint cali- 
bration errors in each filter, and the constraint can arise 
from subsidiary measurements such as on standard stars 
or instrumental calibration (see, e.g., [3l|) that limit the 
zeropoint offsets to lie within AZ G [Z- , Z+] . Again, it 
is the relative zeropoint differences, or colors, that cause 
bias; uniform offsets do not affect cosmology. Because 
the flux from sources at different redshifts peaks in dif- 
ferent wavelength bands, the zeropoint errors will induce 
a redshift dependent error in the magnitude and hence 



a bias in the cosmology parameter estimation. In addi- 
tion, because the use of multiple bands to define the dust 
correction leads to an interdependence of SN at different 
redshifts, a correlated error matrix enters [30| . 

In the two band toy model for dust correction, the 
corrected magnitude m is related to the magnitudes mea- 
sured in two neighboring bands by 



{I + R) rrii — Rm 



i+l 



(16) 



where R is the extinction ratio. We use as the two bands 
the restframe B and V bands for each supernovae, take 
R = 2.1 (somewhat emphasizing the effect), and consider 
only calibration zeropoint error contributions to the dust 
correction, not any intrinsic SN color variation. This 
simple model is sufficient to illustrate the effects of a 
nondiagonal error covariance matrix 



C = BEB^ 

B = {l + R)S,,-Rd.,^,-i, 



(17) 
(18) 



where E is the pre-correction, possibly diagonal, error 
covariance matrix. 

The Fisher matrix is formed using the nondiagonal er- 
ror covariance matrix C and we then calculate the param- 
eter bias due to zeropoint offsets AZk in the magnitudes 
by means of Eq. (|12l) . We consider 8 filters logarithmi- 
cally spaced in wavelength, with centers at Aq (1-l-a*)'''"^, 
for fc = 1 - 8, taking Aq = 4400A and a* = 0.15 so the 
maximum redshift corresponds to 1.66 ^32,] . Our state 
function is AZ, which can be both positive and negative 
within the bounds, and we scan over all possible forms 
within the bounds and analyze the cosmology bias. 

Table |TT] presents the results in the same format as 
the previous case in Table HI However here the steps 
are in band zeropoints not population fractions and the 
locations are listed in terms of the filter numbers. The 
important quantity is AZ^ei between filters (recall that 
an overall zeropoint error has no cosmology effect), and 
needs to be constrained to the ~ 0.01 level. The furthest 
red filters, used for the highest redshift SN, are among 
the most sensitive to bias and should be tightly calibrated 
with instrumental and standard star measurements. 

As in the previous case in Sec. lIVl the shape, i.e. wave- 
length dependence of the zeropoint calibration, is un- 
affected by amplitude of the bounds on the calibration 
state function, but the absolute level is determined by 
the bounds. If we consider twice as large values for AZ, 
then Ax^ just scales with AZ^. 

Any filter zeropoint step affects SN in some redshift 
range, and hence the overall cosmology. Thus an im- 
provement in the knowledge of a filter offset AZk reduces 
the maximum Ax^ bias. Figure [TT] shows the effect of 
more tightly calibrating a given filter^. We see the great- 
est improvement for the end bands - which are used for 



^ We have also carried out PCA on the filter model Fisher matrix 
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Parameters Transitions max Axaz=o.oi A-^rei (Ic) 





2-3, 7-8 


4.07 


0.0098 


Wo 


1-2, 4-5, 7-8 


2.79 


0.012 


Wa 


1-2, 4-5, 7-8 


2.99 


0.012 




3-4, 7-8 


4.12 


0.015 


flm, Wo, Wa 


3-4, 7-8 


4.55 


0.018 



TABLE II: For each set of parameters we consider the form 
of filter zeropoint errors that maximize the bias in terms of 
Ax^ . The transitions column gives the filter transitions of the 
maximizing function in AZ, delivering a maximum bias Ax^ 
scaled to the case where the zeropoint calibration is bounded 
by \AZ\ = 0.01, shown in the next column. The last column 
shows the value of AZrci between two filters that will shift the 
derived cosmology by Icr from the true cosmology. 



the lowest redshift, anchoring SN and the highest red- 
shift, lever-arm SN - and the middle bands, which are 
used near the sensitive z w 1 region. Thus survey design 
that provides particularly comprehensive calibration for 
these bands will see a large payoff in systematics control. 




0.000 0.002 0.004 0.006 0.008 0.010 
AM 

FIG. 11: The maximum Ax^ bias is plotted as a function 
of the bound on individual filter band zeropoint offsets. The 
greatest improvement in systematics control comes from im- 
proving the calibration uncertainty AZ for bands 1, 5, 6, 8 
which correspond to the most cosmologically sensitive redshift 
leverages, near z m 0, 1, 1.7. 



but find little useful information from the modes. The wave- 
length band bin basis is better suited to actual design than say- 
ing, e.g., calibrate the combination of 0.4 times the first filter 
and —0.2 times the second filter etc. 



VI. CONCLUSIONS 

Modeling unknown functions in cosmology is pervasive, 
whether these are functions that carry the physics we are 
directly interested in, such as the dark energy equation 
of state history, or functions describing subsidiary effects 
that we wish to subtract out, such as intermediate as- 
trophysics or modeling of uncertainties. If the functional 
form is assumed, then this becomes parameter fitting or 
"self-calibration" but it is interesting and important to 
investigate the results when any viable function is al- 
lowed. The viability is subject only to some constraints 
placed on the bounds of the function through theoretical 
or measurement input. 

We investigated two main issues for cosmological anal- 
ysis in the presence of unknown but constrained func- 
tions. First, we demonstrated a computationally effi- 
cient, pure and complete method for determining the 
viable space of principal component coefficients, with a 
simple geometric picture in terms of a hyperrectangle in 
the iV-dimensional basis space. This improves on the effi- 
ciency of the previous mode-by-mode method by a factor 
3, while guaranteeing purity, i.e. validity of the selected 
functions. Conversely, compared to the also pure and 
complete direct scanning method, the projected hyper- 
rectangle method gains in efficiency by factors of 10^" or 
greater. 

For many astrophysical problems the orthogonal bin 
basis (in redshift or wavelength) is well suited. We 
evaluated two "real world" systematics issues with this 
method, one in redshift and one in wavelength. The 
first dealt with unrecognized population evolution in a 
subclass of standard candles. Rather than assuming a 
form of evolution we analyzed the effects of all possible 
functional forms lying within some bounds. The results 
provide quantitative guidance to controlling the worst 
cosmology biases arising from the systematic uncertain- 
ties. In particular, we find that the regions z w 0.1 and 
z « 1.0 can benefit most from comprehensive observa- 
tions to limit unrecognized subclasses (at the ^ 0.01 mag 
level); surveys homogeneous over these redshift regions 
have improved control over cosmology misestimation. 

The second application concerned dust extinction cor- 
rections for Type la supernovae, where measurements 
in multiple wavelength bands can fit for dust but also 
correlate supernovae at different redshifts. We analyzed 
the case where systematic uncertainties existed in the 
filter calibrations, bounded by instrumental or standard 
star observations, and propagated all possible functional 
forms into cosmology biases. Again we found which spe- 
cific forms were most damaging and that measurements 
designed to control such errors could remove the worst bi- 
ases. In particular, those bands used for the lowest and 
highest supernovae, and ones relevant around z « 1.0, 
should be most comprehensively calibrated (to the ~ 0.01 
mag level relatively) for a more robust survey. 
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Appendix: Efficient Projection and Boundary 
Definition 

Since the dark energy physics of interest mostly trans- 
lates into descriptions of the EOS in redshift bin space, 
our bounds on the state function are given in this space. 
In this Appendix we will detail, through a geometrical 
understanding, a fast method for analyzing how these 
bounds in bin space may be carried into bounds on co- 
efficients in principal components space. The bounds on 
w{z) directly form the space of the permitted sets of these 
coefficients, but picking a pure and complete set of coef- 
ficients in PC space is not so simple. 

The value for the EOS, say, at a particular redshift z is 
not described by just one mode as in the bin space case, 
but by a linear combination of many modes weighted 
by their respective coefficients. To scan the whole coef- 
ficient space within the bounds given in redshift space 
requires an impractical amount of computational power, 
as discussed in Sec. [TTl We need a more clever method of 
obtaining the desired results. 

For a set of independent state function bounds on {wi}, 
the permitted space in bin space is bounded by a hyper- 
rectangle, or orthotope. The coordinates of the corners 
are given by any combination , , • ■ • , where 
denotes the maximum and minimum values allowed 
for bin i, and N is the number of bins. This fixes the 
2^ corners. This orthotope structure contains all the 
needed information on the boundary of the allowed space 
no matter the basis. 

Let us first address how to determine the viable re- 
gion for any set of principal component coefficients, 
ai, . . . ,am- The permitted space is simply the pro- 
jection of the orthotope onto the subspace spanned by 
these components. The nodes of the boundary in PC 
space are found by projecting the vectors going from 
the origin to the corners of the orthotope. Denoting 
these (A^-dimensional) vectors as Ci and the projected 
(m-dimcnsional) node vectors in the S = [ei, . . . , e^] PC 



space as pi = {pn, . . . 

iPim)} where p^i, . . . ^Pim '^re the 
coordinates of pi with respect to the axes ai, . . . , am, we 
find 

Pi = (S^S)-iS^Ci. (19) 

Consider now the case where we keep only a subset 
of PC modes, i.e. we do not marginalize over the other 
modes but fix their coefficients, e.g. to 0. Information 
is lost by not using the complete basis in the expan- 
sion, but sometimes data or practicalities do not enable 
us to know or measure all modes. Therefore, we must 
sometimes sample with only a subset M of modes. The 
new permitted space of PC coefficients is the intersec- 
tion of the M-dimensional space with the original TV- 
dimensional space. This subspace is completely defined 
by the boundary points, found by looking at the inter- 
sections between the M-dimensional space of interest and 
the A^-dimensional bounded bin space (cf. the black dots 
in Fig. [5]). However, the allowed space is not generally 
a M-rectangle; for example, a plane cutting through a 
cube can have a boundary with six corners unless it is 
specially oriented. This has the consequence that there 
is no A/-dimensional basis that will generally be pure and 
complete (i.e. there are no orthogonal axes spanning the 
space). Only the original bin basis in A^-dimensions can 
provide a pure and complete description of the valid state 
functions. 

To describe the actual procedure for evaluating the al- 
lowed region in PC space we rephrase the issue more 
mathematically. The determination of the coefficient 
bounds is solved by considering the intersection of the 
subspace Ca, spanned by the principal components keep- 
ing the selected number of coefficients constant, with the 
orthotope. In the 3 dimensional case, with a 2 dimen- 
sional subspace, the boundaries are of 3 — 2 dimensions, 
i.e. lines, having parametric form la + (It — la)^ and the 
plane is parameterized by Pq + (Pi — Po)u + (P2 — Po)''^j 
where t, u, v are real numbers, and 1^ and Pj are 
points on the line and in the plane, respectively. Set- 
ting D = [la - Po], A = [I„ - lb. Pi - Po, P2 - Po], and 
L = [t, u, v] we find 

D-^A = L. (20) 

In the higher dimensional case the solution has the iden- 
tical form. Geometrically, I represents vectors going from 
the origin to the corners of the orthotope (so t G [0, 1]) 
and P are non-colinear points in Ca. The matrix manip- 
ulations can be computed easily so solving for the allowed 
region in PC space is highly efficient and quick. 



[1] R. de Putter & E.V. Linder, JCAP 0810, 042 (2008) [2] D. Huterer & G. Starkman, Phys. Rev. Lett. 90, 031301 
iarXiv:0808.0189, (2003) arXiv:astro-ph/0207517, 



13 



D. Huterer & A. Cooray, Phys. Rev. D 71, 023506 (2005) 
arXiv:astro-ph/0404062 

R. de Putter & E.V. Linder, Astropart. Phys. 29, 424 
(2008) arXiv:0710.0373; 

M.J. Mortonson, W. Hu, D. Huterer, Phys. Rev. D 79, 
023004 (2009) arXiv:0810.1744] 

T.D. Kitching & A. Amara, MNRAS 398, 2134 (2009) 
(arXTv:0905.3383 

M.J. Mortonson.' larXiv:09 08.0346 
M.J. Mortonson & W. Hu, ApJ 
[arXiv:0705.TT32] 

E. V. Linder, P hys. Rev. D 79 
larXiv:0812.0370] 
T.D. Kitching, A. Amara, F.B. AbdaUa, B. Joachimi, A 
Refregier, MNRAS 399, 2107 (2009) arXiv:0812.1966 
D. Branch, S. Perlmutter, E. Baron, P. Nugent 
a rXiv:astro-ph/0109070, 
M.Sulhvan 



672, 



023509 



737 (2008) 
(2009) 



et al.. 



ApJ 648, 



868 



arXiv:astro-ph/0605455 
D.A. Howell, M. Sullivan, A. Conley, R.G. 
667, L37 (2007) arXiv;astro-ph/0701912 
T .J. Bro nder et al., A&A 477, 

arXiv:0709 .0859 

D.A. Howell et al., ApJ 691, 661 (2009) arXiv:0810.0031" 
M. Sullivan et al., ApJ 693, L76 (2009) arXiv:0901.2476' 
S. Linden, J-M. Virey, A. Tilquin, A&A 506, 1095 (2009) 



(2006) 
Carlberg, ApJ 
717 (2008) 



[is: 

[19 

[20 
[21 

[22; 

[23; 

[24; 

[25' 
[26' 
[27 
[28' 
[29; 

[30 

[31 

[32; 



[arXiv:0907:4495] 

E.V. Linder, Astropart. Phys. 26, 102 (2006) 
arXiv:astro-ph/0604280 

M.G. Kendall, A. Stuart, J.K. Ord, Advanced Theory of 
Statistics (Oxford U. Press: 1987) 
C. Shapiro, ApJ 696, 775 (2009) |arXiv:0812.0"769l 
S. Dodelson, C. Shapiro, M. White, Phys. Rev.'O 73, 
023009 (2006) arXiv:astro-ph/0508296 
A.G. Kim, E.V. Linder, R. Miquel, N. Mostek, MNRAS 
347, 909 (2004) arXiv:astro-ph/0304509^ 

PASP 118, 205 



T.M. Davis , B.P. Schmidt, Al T 
(2006) , arXiv:astro-ph/0511017 ' 

^J 



Kim, 



663, 



1187 (2007) 



E.Y. Hsiao et al., 

arXiv:astro-ph/0703529 

A. Conley et al., ApJ 664, L13 (2007) (arXiv:0705,03_67J 
M. Kowalski et al., ApJ 686, 749 (2008) " arXiv:0804.4142 
M. Hicken et al., ApJ 700, 1097 (2009) arXiv:0901.4804^ 
S. Nobih et al., ApJ 700, 1415 (2009) arXiv:0906.4318 ' 
J. Nordin, A. Goobar, J. Jonsson, JCAP 0802, 008 (2008) 

arXiv:0801.2482 

A.G. Kim & R. Miquel, Astropart. Phys. 24, 451 (2006) 
arXiv:astro-ph/05082 52 

C.W. Stubbs et al., A SP Conf. Series 364, 373 (2007) 
arXiv:astro-ph/0609260 
G. Aldering et al., arXiv:astro-ph/0405232| 



